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Abstract 

In this paper, we investigate the effect of spatial dispersion on a double-lattice metamaterial 
which has both magnetic and electric response. A numerical scheme based on a dipolar model is 
developed to extract the non-local effective medium of the metamaterial. We found a structure- 
induced bianisotropy near the resonance gap even if the artificial particles are free from cross- 
coupling. A cross-coupled resonance also results from spatial dispersion and it can be understood 
using a Lorentz frequency dispersion model. The strength of the cross-coupled resonance depends 
on the microstructure and it is larger at a higher normalized frequency. 

PACS numbers: 78.20.Ci,73.20.Mf,42.70.Qs,78.67.Bf 
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I. INTRODUCTION 



Metamaterials consisting arrays of artificial particles usually work in the frequency regime 
in which the free-space wavelength is larger than the inter-particle distance by a factor 
around 10. Their electromagnetic properties can therefore be understood by treating a 
metamaterial as a homogeneous medium. More often, a local effective medium of permit- 
tivity (e e //) an d permeability (/i e //) is assumed to be valid. There are several approaches in 
obtaining the local effective medium of a metamaterial. One approach is the S-parameter 
retrieval method.- It obtains the material parameters by comparing the transmission and 
the reflection coefficients of a slab of metamaterial and the corresponding homogeneous slab. 
Another approach relies on a homogenization theory which averages the microscopic field by 
line and area integrals.-^ Both approaches can cater for complex geometry and microstruc- 
ture of the metamaterial and they work very well at frequencies far away from the resonance 
gap where spatial dispersion becomes important so that higher order eigenmodes (including 
the non-propagating ones) can be excited inside the metamaterial.- In this work, we will 
take another approach to deal with the spatial dispersion at a frequency near the resonance 
gap. Instead of picking or assuming a dominant eigenmode within the metamaterial, we 
will find the non-local effective medium directly for an arbitrary Bloch wave vector q. For 
instance, the permittivity e(u, q) now depends on both the angular frequency uj and the 
wave vector. 

Apart from working near the resonance gap, some metamaterials like the wire medium 
the three-dimensional wire mesh 8 and the structured metal surface^ also indicate that 
spatial dispersion (the non-local effective medium) can become important in order to have 
an accurate description of the metamaterial. In this paper, we will establish a numerical 
scheme in getting the non-local effective medium and some effects of spatial dispersion on 
metamaterial will be examined. 

Section [Til outlines the numerical scheme in obtaining the non-local effective medium. We 
will employ a dipolar model for simplicity. It has the advantage that analytic formulas can 
be obtained in some cases. Moreover, all the material parameters, two tensors {e and //) 
together with two pseudo tensors (£ and () can be obtained easily in the whole frequency 



regime. All the four constitutive tensors will be extracted as Ref. 



10 and 
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pointed out that 



a metamaterial should be described as a bianisotropic medium if the artificial atoms suffer 
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from cross-coupling (electric field generates magnetic dipole and/or magnetic field generates 
electric dipole). 

In Section 1111} the local effective medium of a double-lattice metamaterial will be first 
examined. One sublattice consists of electric artificial atoms and another sublattice consists 
of magnetic artificial atoms. The condition that the local effective medium becoming free 
of bianisotropy will be discussed. 

In the last section, the non-local effective medium of the double-lattice metamaterial will 
be examined. The relationship between the non-local effective medium to the corresponding 
local one will be discussed. On the other hand, a Lorentz-type model is suggested for fitting 
the frequency dispersion of the non-local effective medium. Some effects of the spatial 
dispersion, including the structure-induced bianisotropy and the cross-coupled resonance, 
will also be discussed. 

II. FORMULATION 

Here, we give an outline of the numerical scheme in obtaining the non-local effective medium 
using a dipolar model. The artificial particles ( split-rings, short wires, or electric atoms,— 
etc.) are placed in an infinite lattice of lattice vectors R in vacuum (of wavenumber k , wave 
speed c , and intrinsic impedance rj ). In each unit cell labeled by R, there are n particles 
(indexed by j = 1, n) at positions R + Tj. In the dipolar limit, each artificial particle has 
an electric dipole moment p (R + Tj) and a magnetic dipole moment m (R + Tj) and they 
are governed by the following equation (with time dependence factor e~ lu)t and S. I. units): 

C p (R + Tj) 

_ m(R + Tj) 

where Ej Qca j (R + Tj) and Hj ocg j (R + Tj) are the local electric and local magnetic fields at 
the center of the j-th particle. These local fields include the fields radiated by the external 
sources and the scattered fields from all the other particles. The 6-by-6 matrix ctj represents 
the polarizability of the j-th particle. 

In this model, every artificial particle is replaced by an equivalent dipolar point particle. 
In order to extract the non-local effective medium of this crystal at an arbitrary wave vector 
q, we drive the crystal by the external fields in the form 



E local ( R + r i)A?o 
H local ( R + T i) 



3 



(2) 



E ext ( r ) = E ext ex P (*q ■ r ) > 

H ext ( r ) = H ext ex P ( ic l • r ) • 

These external fields are in fact generated by an electric and a magnetic polarization of the 
same spatial dependence. The crystal responds to it by generating dipole moments in the 
form 



(3) 
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which satisfies the multiple-scattering equation: 
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H ext 



(4) 



with V being the volume of a single primitive unit cell. The 6n x 6n transition matrix T is 
given by its inverse: 



(T- 1 ) . ., (q) = VajHtf - G jr (q) , 
where the lattice Green's function Gjf is defined by 



(5) 



G 
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r=R+Tj-T,,^0 



e- iq r G (r). 



(6) 



The 6-by-6 matrix Go (r) which relates a dipole moment to its radiation field, together with 
its Fourier transform Go (q) are given explicitly in Appendix [A] The lattice Green's function 
can be evaluated by the Ewald sum technique. For mathematical convenience, we also define 



ok 3 V 

G a> (q) + -i^Ui? = Go (q) + 



[,r (q) ijf (q) 
-7«'(q) 3ii'(q) 



(7) 



where 7 is the 6x6 Identity Matrix. The tensors /?,- •/ (q) and 7,y (q) are periodic functions 
in the reciprocal space. They can be characterized according to the point group^ associated 
with the lattice structure at a particular wave vector q together with the requirement that 
the matrix transpose (fixed j and f) of 7jy (q) satisfies 
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ijj> (q) T = -ijj> (q) , ( 8 ) 

which is already implied from Eq. ([6]) with Eq. (IA2I) . Table [I] lists the forms of the tensors 
(3jj, (q) and 7^, (q) for some common single lattices and double lattices. The matrices 
are listed with respect to the Cartesian basis vectors {x, y, 2} where z is defined to be the 
direction along q in our convention. 
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TABLE I: The tensor forms of the lattice Green's function for different single and double lattice 
structures along different directions of the wave vector. 

In the next step, the dipole moment distributions numerically solved from Eq. (J3J) are 
then macroscopically averaged. The procedure we choose here is to apply a low pass filter 
such that all the Fourier components outside the first Brillouin zone are eliminated. In 
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essence, it means that all the microscopic fields are ensemble averaged by translating the 
whole crystal with an arbitrary distance while the external fields remain unchanged. Thus, 
the macroscopic polarizations are written as 



c (P) (r) = c (P) e iqr = I 

1 n 

(M) (r) = (M) = - £ 

i=i 

The macroscopic polarizations can now be formally written as 



V qr , 



m,-e-' q - 7 V q - r . 



co (P) 


= CO (q) 


Eext/^o 








H ext 



where the averaged transition matrix is given by 



(9) 



(10) 



(11) 



CO(q)= £7>(q)- 

Eq. fjiup summarizes the response of the crystal to the external fields. Then, the macroscopic 
fields are constructed from the sum of the external fields and the secondary radiation by 



(12) 
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Finally, by combining Eq. (jTOjl and Eq. ( 1121) . we obtain the constitutive relationship 
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e(q)-7 ?(q) 




(E) /770 


_ (M) 
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(13) 



where the four constitutive tensors are given by 



(T)- 1 (q) + G (q). 



(14) 



e(q)-/ C(q) 

C(q) M(q)-? 

In general, apart from the permittivity tensor e (q) and the permeability tensor // (q) , we 
also obtain the off-diagonal tensors £ (q) and £ (q) which represent the bianisotropy of the 
metamaterial. We will see in the next section that the bianisotopy should not be neglected 
in considering the non-local effective medium. All the four tensors are functions of the 
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frequency ui and the wave vector q. For convenience, we will omit writing u inside the 
parenthesis for the constitutive tensors. Moreover, we will further omit the dependence of 
q if the local limit q = is taken, i.e. t (q = 0) is abbreviated as t. 

We note that in our scheme of defining the non-local effective medium at a particular 
q, a fixed form of external fields (Eq. (121)) is used. There are alternative methods which 
directly deal with the macroscopic fields without considering the external fields.— 1 ^ 

III. THE LOCAL EFFECTIVE MEDIUM - DOUBLE LATTICE 

In this work, we will concentrate on the double lattice structure which is a common type 
of microstructures for metamaterials in obtaining a negative refractive index. The double 
lattice consists of two types of artificial particles. One type of the particles (labeled by j = 1) 
has a resonating electric response (e.g. a short wire) and another (labeled by j = 2) has a 
resonating magnetic response (e.g. a split-ring). In particular, we would like to study the 
importance of the electro-magnetic coupling between the electric and the magnetic particles 
when we have a crystal structure having both electric and magnetic response. It is a general 
feature of metamaterials. 

First, we would like to restrict our discussion to the class of metamaterials whose local 
effective medium is free from bianisotropy. It can be achieved through a careful design of 
the artificial particles such that the particles are free from cross-coupling.— 1 ^ Moreover, the 
lattice structure has to be carefully chosen as well. In fact, all the double lattice structures 
listed in Tabled] have vanishing 7--, at q = 0. From Eq. (jSJ), Eq. (TTTT) and Eq. (}Tlj) . it 
means that 



For simplicity, we further assume that the electric particle can only be electrically polarized 
along the x-direction with polarizability a e (c pi x = a e Ei oca ^ x {r-i) / 'rj ) , the magnetic particle 
can only be magnetically polarized along the y-direction with polarizability a m (vn^y = 
a m H 'i oca i ,y(r2)) and all the polarizabilities along the other directions are assumed to be 
negligible. In this case, the relevant basis vectors can be reduced to {c pi x , m 2y } in Eq. (j4j). 
Eq. (fTUl) in the local limit can then be written as 




(15) 
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is now written in the {cq (P x ) , (My)} basis. The macroscopic electric polarization (P) is 
only contributed from the first particle while the macroscopic magnetic polarization (M) is 
only contributed from the second particle. Therefore, the local effective medium is governed 
by 



1 



e T 



— + fl T = Va; 1 + iklV/ (6tt) 
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(18) 



+ f3 T = Va^ + ik d V/(Qir) 



fMr-1 

Because of the carefully chosen microstructure, the local effective medium has exactly the 
same form to a single lattice of one type of particles.— This is the ideal case that the electric 
and the magnetic particles cannot "see" each other. Here, (3t is the local limit (q = 0) 
of j3r (q) defined in Eq. (JTj). It is a function of the normalized frequency uia/ (2nco) for 
the lattice structure and it can be well approximated by a polynomial at small frequencies. 
Table [IT] gives the polynomial expansion of (3t for various lattice structures. Note that 
for cubic lattices, Pt 1/3 at the long wavelength limit so that Eq. (fT8j) returns to the 
Clausius-Mossotti equation.—^ 





0T 


SC 


1/3 - 5.970 2 + 11.8ft 4 


FCC 


1/3 - 2.40ft 2 + 1.72ft 4 


BCC 


1/3 - 3.82ft 2 + 4.38ft 4 



TABLE II: Polynomial expansion of (5t for different single lattices. The normalized frequency is 
ft = u>a/ (27rco) and a is the lattice constant. 



One popular form of the frequency dispersions of the local effective medium is given by 
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Here, Eq. f|T9l) is still called the Lorentz-model dispersion although the resonance term in 
permeability is modified so that these frequency dispersion forms give us the correct low 
frequency behavior (including the frequency regime near resonance) when there are only 
one electric resonance and one magnetic resonance.— One point to note here is that for a 
passive particle, we have to satisfy Im (l/a e ) < — k$/ (Qir) for the electric polarizability and 
also a similar inequality for the magnetic polarizability— It guarantees a positive imaginary 
part for both the local permittivity and the local permeability according to Eq. (fl8l) . i.e. T e 
or r, m should be zero for a lossless medium or positive for a lossy medium under the dipolar 
model. 

In the following, we would like to specify the local effective medium {er, [it} directly 
instead of specifying the polarizabilities {a e , a m }. As an example, we consider the case 
that uj (2tt) = 8GHz, A = 0.89, uj m / (2tt) = 8.5GHz, B = 0.128 and T e = T m = 0. The 
underlying microstructure is assumed to be a CsCl structure of lattice constant a = 5mm. 
The polarizabilities can be found by Eq. (1181) whenever a microscopic description is neces- 
sary. Fig. 1(a) and (b) show the local effective permittivity and the local permeability. Fig. 
1(c) shows the crystal being viewed from the direction [001] or [010]. The uniaxial electric 
particle (green color) can be polarized along the [100] direction and the uniaxial magnetic 
particle can be polarized along the [010] direction. We will start from this configuration in 
the next section when we consider the non-local effective medium. 

It is more convenient to specify the metamaterial by its local effective medium directly. 
However, a note has to be taken for the freedom of the parameters in the frequency disper- 
sion. In particular, for a non-absorptive particle, the polarizability can be written in terms 
of the dipolar scattering phase shift (a real number) by 



where a = e,m for "e" wave scattering (non-zero local electric field with vanishing local 
magnetic field at the particle) or "m" wave scattering (non-zero local magnetic field with 




(20) 
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FIG. 1: (Color online)A double-lattice configuration specified by Eq. (|19p with u e / (2ir) = 8GHz, 
A = 0.89, uj m j (2ir) = 8.5GHz, B = 0.128, T e = T m = and lattice constant a= 5 mm. (a) local 
permittivity; (b) local permeability; (c) The underlying CsCl structure viewed along the [001] and 
the [010] direction. The uniaxial electric particle (green color) can be polarized along the [100] 
direction. The uniaxial magnetic particle (blue color) can be polarized along the [010] direction. 

vanishing local electric field at the particle). Due to the finite volume of one single unit 
cell, the phase shift d[ a ^ cannot vary arbitrarily against frequency. Suppose a sphere of 
radius r s can completely enclose one artificial electric particle, the restriction on 8^ can 
be written in terms of the time-averaged total energy <§ within the enclosing sphere (for a 
quasi-monochromatic electromagnetic field) as: 



(21) 



jr( e ) i , . 

^<g oc^^l in + 29 2 - 29 4 - cos 2 ( 9 + 6[ e) ) 

duj 29 A V / 

-29 sin 2 (W<^ e) )) > 0, 

where 9 = ur s /c . See Appendix |B] for the derivation. This inequality must be satisfied. 
Otherwise, we can extract energy from the artificial particle without first pumping energy 
to it. 

For our example, the polarizabilities are obtained through Eq. (1181) and the time-averaged 
total energy for "e" wave scattering (in an arbitrary unit) is obtained from Eq. fl2"01 and 
Eq. (12~T1) . Fig. [2] shows the results for two different cases. The solid line shows the case 
with r s = 2mm. The total energy within the sphere is positive in the whole frequency 
regime. The dashed line shows the case with a smaller radius r s = 1.4mm. In this case, the 
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FIG. 2: The time-averaged total electromagnetic energy within a sphere of (a)r s = 1.4mm and (b) 
r s = 2.0mm enclosing the electric particle for "e" wave scattering. 

energy becomes negative which is not valid. Therefore, in order to design our metamaterial, 
the electric artificial particle must have a certain minimum size. On the other hand, if we 
already know the size of the electric artificial particle, it imposes a maximum on the electric 
resonance strength A (in Eq. (j!9p ) that the metamaterial can have. 

IV. THE NON-LOCAL EFFECTIVE MEDIUM DOUBLE LATTICE 

In this section, we investigate how the effective medium changes when the wave vector 
is deviated from the Brillouin zone center. In fact, the tensor 712 (q) does not vanish in 
general. For the double lattice we considered in the last section with q along direction [111] 
or [100]/[010]/[001], Eq. (HOD can still be written in the form of Eq. (HH where 



This 2-by-2 matrix is in the {co (P x ) , (M y )} basis with z being the direction along q, x 
being the direction the uniaxial electric particle can be polarized and y being the direction 
the uniaxial magnetic particle can be polarized. Therefore, by comparing to Eq. (|T7|) . the 
non-local effective medium can be written with respect to the local effective medium by 



((T^ + Go) (q) 




(22) 



11 



-1 -1 



er(q)-l -£ (q) 

-C (q) (q) - i 



^ - A/3 T (w, q) 712 (w, q) 

7i* 2 (<*>, q) - (w, q) 



(23) 



where er(q), /ij-(q), £ (q) and £ (q) are the elements of the four constitutive tensors in 
y, z} basis: 
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(24) 



and A(3t (k>, q) is defined by A/3r (o>, q) = (3t (o>, q)— /3y (w, q = 0). e (q) in the y/z direction 
and Ju (q) in the x/z direction have the value one since we have neglected the polarizabilities 
in these directions for mathematical simplicity. 

Now, we use Eq. ( )23l) to find the non-local effective medium for the previous example 
shown in Fig. 1. In this case, we plot e*r (q), fPr (q) an d £ (q) along the T to Z[001] direction 
in Fig. [21 The profile of (q) and fiT (q) are shown in Fig. [3] (a) and (b). For every wave 
vector q, we can see that €t (q) and \xt (q) have the similar frequency dispersion as the 
local medium correspondence and fiT, only with the resonating frequency being shifted. 
This qualitative behavior is expected from Eq. (123]) by neglecting the presence of 712 (q). 
However, with careful examination, we can still see that when q deviates from the T point, 
tT (q) also diverges at the magnetic resonating frequency. Moreover, we have non-vanishing 
£ (q) near the electric or the magnetic resonating frequency. These result from the presence 
of the term 712 (q) which means the cross-coupling between the electric and the magnetic 
field due to spatial dispersion. 

For the current configuration, the effect of 712 (q) is rather weak. However, it can become 
more prominent in other cases. As an example, for the same CsCl lattice structure, we 
now consider the wavevector q along the T to R[lll] direction and the particles are now 
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FIG. 3: The non-local effective medium from the V to the Z[001] point for the configuration shown 
in the caption of Fig. 1: (a) ej- (q); (b) fir (q) and (c) £ (q). White color denotes positive and 
black color denotes negative values. 

orientated so that the uniaxial electric particle can be polarized along the [-112] direction 
while the uniaxial magnetic particle can be polarized along the [1-10] direction. This second 
configuration is shown in Fig. HI ey(q), f^r (q) along the T to R[lll] direction are now 
plotted in Fig. [5]. In Fig. H^a) and (b), we also plot ex (q) and \it (q) against frequency 
at exactly the R point. From the results, the resonating frequencies depend on q and both 
€t (q) and \xt (q) n o w clearly diverge at two separate frequencies (instead of one) for every 
non-zero q due to the appearance of the term 712 (q). The two resonating frequencies U\ (q) 
(the lower one) and 0J2 (q) (the upper one) can be obtained from Eq. ( 1231) and are governed 
by 

( — T " A/5t (q)) ( l —r - Af3 T (q)) 

\e T -l J \ht-1 J (25) 

- 1712 (q)| 2 = 0. 

If the magnetic resonance and the electric resonance are far apart in frequency or the term 
7i2 (q) can be neglected, it decouples to 
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FIG. 4: (Color online)A double-lattice configuration specified by Eq. (fl~9|) with uj e / (2ir) = 8GHz, 
A = 0.89, Lo m / (2tt) = 8.5GHz, B = 0.128, T e = T m = and lattice constant a= 5 mm. The 
underlying CsCl structure is viewed along the (a) [111] and (b) [0-10] direction. The uniaxial 
electric particle (green color) can be polarized along the [-112] direction. The uniaxial magnetic 
particle (blue color) can be polarized along the [1-10] direction. 
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FIG. 5: The non-local effective medium from the T to the R point for the configuration shown in 
the caption of Fig. HJ (a) et (q); (b) fix (q) and (c) £ (q). White color denotes positive and black 
color denotes negative values. 
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FIG. 6: The non-local effective medium at the R point for the configuration shown in the caption 
of Fig. m (a) e T (q); (b) /x r (q) and (c) £ (q). 



- A#r (q) = 0, 

ey — 1 

which determines the electric resonating frequency and 
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AP T (q) = 0, 



(27) 



(It - 1 

which determines the magnetic resonating frequency. It gives us two separate "bands" of 
resonating frequencies. The term A(3t (q) gives us the q dependence of the electric or the 
magnetic resonating frequency (shown in the previous example in Fig. [3]). However, in the 
current example, these two bands hybridize with each other so that e-r (q) or /Mr (q) has two 
resonating frequencies instead of one. It is called the cross-coupled resonance here. 

Due to the same reason, the 712 (q) term also causes a structure-induced bianisotropy 
when we bring the electric and magnetic resonances near to each other in frequency. This 
bianisotropy is induced by the structure instead of the cross-coupling effect of the artificial 
atoms. The £ (q) from the T to the R point is plotted in Fig. [5]^c). It also diverges at U\ (q) 
and u>2 (q). 

We have used a Lorentz- model dispersion (Eq. ( |T9i) ) for the local effective medium. In 
fact, er(q), fJ-r (q) and £ (q) can be numerically fitted very well by extending the same 
model: 
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In the current example at the R point, the expressions obtained by setting u\ (q) / (2tt) = 
8.57GHz, uj 2 (q) / (2tt) = Q.llGHz, A x (q) = 0.13, A 2 (q) = 0.55, B 1 (q) = 0.10, B 2 (q) = 
0.03, C (q) 27r = O.lbGHz^ 1 and Ti (q) = T 2 (q) = can fit the actual results very well. 
The fitted result has no noticeable differences from the results shown in Fig. El 

The extended Lorentz-model dispersion is valid even for an absorptive system. Suppose, 
now we add some absorption to the system by setting Y e / (2tt) = 0.05GHz and T m / (2tt) = 
0.02GHz in Eq. f|T9|) . the resultant non-local effective medium at the R point is shown 
in Fig. [7] . The non-local effective medium can still be fitted by the extended Lorentz 
model but with non-zero Ti (q) and Y 2 (q). In this case, we set Ti (q) / (2n) = 0.025GHz 
and T 2 (q) / (2tt) = 0.0A5GHz while all other parameters remain the same. The fitted 
expressions also show no noticeable differences to the results shown in Fig. [7] . 

From Fig. [3, we see that the imaginary part of both (q) and fir (q) are positive. For 
a bianisotropic medium, there are no general requirements that the imaginary part of the 
permittivity and the permeability should be positive. However, within the dipolar model, 
from Eq. (1231) with j3r (uj, q) being real, it is easy to show that both (q) and fiT (q) 
have positive imaginary part if the local ex and \xt have positive imaginary part, i.e. Ti (q) 
and T 2 (q) should be positive. In our approach, we parameterize the effective medium by 
using q such that q is always a real vector within the first Brillouin zone. If, on the other 
hand, q is allowed to be complex (e.g. traveling along the dispersion of an eigenmode 
and becoming complex within the resonance gap where no propagating eigenmodes can be 
found), a negative imaginary part in er (q) or \lt (q) may appear around a cross-coupled 
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resonance.— 

The structure-induced bianisotropy £ (q) appears even if the local effective medium does 
not suffer bianisotropy. £ (q) is generally present and it is purely real for a non- absorptive 
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FIG. 7: (Color online)The non-local effective medium at the R point for the configuration shown 
in the caption of Fig. Hwith T e = 0.05 and T m = 0.02: (a) e T (q) ; (b) /i T (q); (c) £ (q). The blue 
solid line shows the real part and the gray dashed line shows the imaginary part. 

system in our example. It disappears in the local effective medium limit only if the particles 
do not suffer cross- coupling and the particles are placed in a carefully chosen lattice. For 
example, if the sublattices of the two types of particles are displaced in the z direction 
(along q) with respect to each other, the 2D symmetry perpendicular to q does not change 
so that the constitutive tensors can still be described by Eq. ([21]). However, £ (q = 0) will 
not vanish and it is a purely imaginary number. On the other hand, if the lattice is not 
displaced but the particles suffer cross coupling, £ (q = 0) also does not vanish and it is 
purely imaginary.— 

The structure-induced bianisotropy and cross-coupled resonance are both due to the term 
7i2(q). In fact, with a smaller lattice constant, the interaction becomes smaller. Fig. [8] 
shows the €t (q), fix (q) and £ (q) along the T to R[lll] direction when the lattice constant 
a is changed a smaller value 2.5 mm. The splitting between the two resonating frequencies 
is smaller than in Fig. [5j On the other hand, if the lattice constant becomes larger, the term 
7x2 (q) cannot be neglected in finding the effective medium or in finding the dispersion. 

We have investigated the constitutive tensors of the metamaterial. After finding the 
effective medium, the dispersion diagram can also be found. The dispersion relation of the 
effective medium is governed by 
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FIG. 8: The non-local effective medium from the T to the R point for the configuration shown in 
the caption of Fig. d]with lattice constant a changed to 2.5 mm: (a) er(q); (b) fix (q) and (c) 
£ (q). White color denotes positive and black color denotes negative values. 




FIG. 9: Dispersion diagram (open circles) for the configuration shown in (a) Fig. 1 along the TZ 
direction and the configuration shown in (b) Fig. H] along the TR direction. The black dashed line 
represents the corresponding dispersion obtained from the local effective medium. The blue solid 
line represents the perturbed dispersion from the local one. 
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det((T)- 1 (q))=0. 

After some algebra, we have 



(31) 



det ((T- 1 ) (q)) 



x 



q2 ~ kl C\2\ 
|g + fcoe(q)| 2 -fco 2 e r (q)/i r (q) ^ J 

(er(q)-l) 0*r (q) - 1) - If (q)| 2 
and the dispersion relation is therefore 

|g + ^(q)| 2 = A: 2 e r (q)/iT(q). (33) 
There is an extra term (q) in the dispersion relation comparing to the one obtained from 
the local effective medium: q 2 = k^x^T- Unlike the case of the local effective medium which 
serves as an approximation of the actual band structure, the dispersion relation shown in 
Eq. fl33|) can be used to find the band structure which shows no noticeable differences to 
the actual band structure if we solve it directly without first finding the non-local effective 
medium. For band structure of a single lattice consisting artificial particles of only electric 



or only magnetic response, see Ref. 
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171 and Ref. 

The dispersion diagram along the direction TZ for the first configuration (Fig. 1) and the 
dispersion diagram along the direction TR for the second configuration (Fig. H|) are shown 
in Fig. [9(a) and (b) respectively. They are plotted in open black circles while the dispersion 
diagram obtained from the local effective medium is plotted in dashed line. Along the TZ 
direction for the first configuration, they look very similar especially near the Brillouin zone 
center. In fact, the formation of the non-local bands can be understood from a perturbation 
point of view by assuming small A/3t (q) and 712 (q). For a fixed q, we can approximate 
the change in square frequency from the local to the non-local dispersion by 

d -%^Akt = (e T - 1) fa - 1) - kl) 

OK o (34) 

xdet((T)- 1 (q)), 

which is evaluated along one of the bands obtained from the local effective medium. (The 
proof is not shown here.) By substituting Eq. (123]) into it and neglect the second order 
effect of Pt (q) an d 7 (q), we have 
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+ Re (712 (q)) 2 v /e T /iT (e T ~ 1) (/*r - 1) • 

This perturbed dispersion is plotted in Fig. [9] in solid line. We can see that it approximates 
the actual non-local dispersion very well in the TZ direction of the first configuration. In 
the TR direction, the first band (having positive refractive index in the local limit) and the 
second band (having negative refractive index in the local limit) hybridize with each other 
to form the actual band. 

For every mode in the dispersion diagram, we can characterize it by defining the refractive 
index n of the mode using q 2 = n 2 kQ and defining its characteristic impedance Z to be the 
ratio between the tangential E and the tangential H field. They are governed by 



' + Vl' (36) 
nZ _ (q) 

1+f (q) ko/q' 

For example, the n and Z along the band in the TZ direction in Fig. [9] are plotted as solid 
lines in Fig. [101 for the double- negative band. The corresponding values for the local effective 
medium (n 2 = 6t^t,Z = a/ ^t/^t) are also shown in dashed lines in the same figure. We 
can see that the local effective medium works very well for the propagating bands unless it 
is near the Brillouin zone edge and near the resonance gap. 

We have examined the non-local effective medium governing the transverse waves. In 
fact, if the artificial particles are oriented so that they can only be polarized in a single 
direction along the z-axis, the metamaterial supports longitudinal waves whose propagating 
direction is also along the z-axis. We can do the similar analysis to obtain the longitudinal 
permittivity £i(q) and the longitudinal permeability /ir,(q). In this case, the electric and 
magnetic response are decoupled, e^(q) and /iz,(q) can also be fitted by using a Lorentz- 
model dispersion: 

e L (q) Rs 1 , 37 

uj 1 - uj 2 (q) + iujT e (q) 

and 
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FIG. 10: Solid line: characteristic impedance and the refractive index of the double-negative band 
in the TZ direction. Dashed line: corresponding values for the local medium. 

^ (q) « 1 - -2 — Pr^r i v ( 38 ) 
w - (q) + «^r m (q) 

The electric longitudinal mode (or the electric bulk plasmon with macroscopic E field along 
q) has a dispersion relation governed by 

e L (q) = 0, (39) 

and the dispersion relation for the the magnetic longitudinal mode (or the magnetic bulk 
plasmon with macroscopic H field along q) is governed by 

Mq) = o. (40) 

With spatial dispersion such that ex,(q) and /xz,(q) depends on q, the band for either the 
electric or magnetic longitudinal mode becomes dispersive to form a narrow band instead of 
a flat line in the local medium picture. 

V. CONCLUSION 



We have established a numerical method based on a dipolar model to obtain the non-local 
effective medium of a metamaterial. In particular, we have concentrated on the double lattice 
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structure in which one sublattice holds the electric artificial atoms and another sublattice 
holds the magnetic artificial atoms. Once the dipolar scattering properties of the atoms 
or the local effective medium parameters are specified, the non-local effective medium with 
all the four constitutive tensors can be obtained on the whole frequency regime including 
frequencies near the resonance gap. We found that the metamaterial should be regarded as 
bianisotropic near the resonance gap even if the artificial atoms do not suffer cross-coupling. 
In this case, the cross-coupling comes from spatial dispersion and it also induces a cross- 
coupled resonance for metamaterials having both electric and magnetic resonances. Within 
our model, the cross-coupled resonance can still be understood using a Lorentz dispersion 
model. The effect of the cross-coupling from spatial dispersion depends on the microstructure 
and it is larger at a higher normalized frequency (or a larger lattice constant). 
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APPENDIX A: THE GREEN'S TENSOR 

The radiation field from an electric dipole p together with a magnetic dipole m at the 
origin in the vacuum is 




G (r) • 



c p 



(Al) 



H(r) 



m 



where the G matrix is defined by 





ik X7 x G (r) 




(A2) 



with the vacuum Green's tensor Go defined by 




(A3) 



where I is the 3-by-3 identity tensor. 

The Fourier Transform of G (r) is given by 
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Go (q) = J G (r) exp (-iq • r) d 3 r 



k^Go (q) -fc q x G (q) 
k q x G (q) k 2 Q G (q) 



(A4) 



where 



^° ^) = ^Ta| " <?<?) " 99- ( A5 ) 

APPENDIX B: RESTRICTION ON THE SCATTERING PHASE SHIFT OF A 
SINGLE PARTICLE 

Here, we will consider the restriction on the polarizability for a lossless classical uniaxial 
particle which can be described by its permittivity distribution e s (r) and its permeability 
distribution /i s (r). The particle has a finite radius r s such that it is vacuum outside the 
sphere of radius r s centered at the origin. Without losing generality, we consider the situation 
that "m" wave (with angular momentum I = 1 and m = 0) is scattered from the particle 
which can be magnetically polarized along the z-direction but not the other directions. We 
also assume that the particle has no cross-coupling, i.e. the electric dipole moment is zero 
when the local electric field is zero here. The total wave outside the particle can be written 
in the following form: 

E=^rxVF Hra= ,(f) (Bl) 

where 

= ! (K (0) + exp (2^! m) ) h (9)) , (B2) 

with hi (6) being the out-going spherical Hankel function, and Y\ m (r) being the spherical 
harmonic. b\ m \ the phase shift in scattering, is a real number for a lossless particle and it 
is related to the magnetic polarizability through 

ik ° a ™ - ew ( 2 i5[ m) ) - 1. (B3) 



3tt 
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On the other hand, from the Maxwell equations, we can prove the following identity 



1 f dujejr) a dujfxjr) 2 3 _ 
AJ duj duj 

— 

The volume and the surface integrals are taken to be the volume and the surface of the 
sphere of radius r s which encloses the particle completely. By substituting Eq. (IB lft into 
Eq. (1B41) . we obtain 



I/ £o ^|E| 2 + , ^|H|VV = 

4 7 a- a. 



4u; 2 /z V ^ ^ 

where = ur s /c. Finally, we recognize the l.h.s. of Eq. (IB5P is just the time-averaged 
total electromagnetic energy ($) within the sphere for a narrow band signal centered at 
the angular frequency uj. By substituting Eq. (1B2I) into Eq. (1B5I) . we can express the 
time-averaged total energy in terms of b~\ ■ 



^ L(i + 2^ 2 - 2£ 4 - cos 2 f + 5™) 

dio 20 A V / 



dS^ 1 

-I I ■ . ■ ^17 — 1\1 — COS i 17 - () , | 

(B6) 

-29 sin 2 (0 + <^ m) J), 

which must be a positive number. The derivation for the electric scattering phase shift is 
similar. Moreover, the derivation also applies for an isotropic particle instead of an uniaxial 
particle being discussed here. 
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